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Abstract 

We  determine  the  phase  separation  and  the  spinodal  of  binary  liquids  with 
anisotropic  interaction  between  the  molecules.  The  results  show  the  influence  of  the 
anisotropy  parameters  on  the  position  of  the  spinodal  curve.  The  model  is 
interpreted  as  a  decorated  lattice  transformation  of  the  Be'.he-Peierls  solution. 
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Since  supercooling  is  most  easily  obtained  in  glasses  and  organic  liquids 
such  as  alcohols,  we  try  to  establish  the  influence  of  anisotropic  interaction 
between  molecules.  One  of  the  simplest  ways  to  introduce  anisotropy  is  to 
assume  that  molecules  have  a  finite  number  of  contact  points,  some  of  these  are 
strongly  adhering  to  similar  points  on  other  molecules  and  other  points  are  not 
strongly  adhering.  The  last  points  may  even  have  a  repulsive  interaction.  In 
order  to  give  meaning  to  the  model  we  assume  that  special  adherence  is 
associated  with  II-bonding  ,  non-special  adherence  with  van  der  Waals  forces  and 
repulsion  may  be  due  to  protruding  parts  of  the  molecule. 

This  model  was  proposed  some  time  ago  by  Barker  and  Fock  to  explain  the 
lower  critical  point  that  occurs  in  the  phase  separation  in  binary  liquids. 

In  this  note  we  point  out  that  about  the  same  solution  of  the  problem  as 
was  obtained  by  Barker  and  Fock,  can  be  found  by  a  different  method.  The 
method  enables  us  to  map  the  solution  on  a  standard  Bethe-Peierls  solution. 

On  the  basis  of  this  we  compute  the  spinodal  and  find  that  anisotropy 
seems  to  "retract"  this  curve  closer  to  the  center.  If  the  number  of  repulsive 
contact  points  dominates  the  spinodal  will  behave  similar  to  a  isotropic  binary 
phase  separation.  Barker  and  Fock1  assume  that  a  molecule  has  two  kinds  of 
contact  points.  For  the  time  being  we  will  assume  that  there  is  one  contact 
point  of  one  kind  (hereafter  called  the  "special"  contact  point)  and  all  the 
other  contact  points  of  the  other  type  (the  "normal"  ones). 

The  interaction  energies  between  two  different  types  of  molecules  will 
depend  cm  which  of  the  contact  points  is  involved  in  forming  the  "bond". 
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Let  us  denote  the  two  types  of  molecules  by  A  and  B  each  of  which  can  be 
oriented  in  n  directions  (absolute  directions  in  space).  This  can  conveniently  be 
described  by  an  arrow,  where  the  tip  of  the  arrow  denotes  the  special  contact  point. 


Thus  for  A  molecules: 
B  molecules: 


r  / 

f  / 


Barker  and  Fock  took  for  n  the  coordination  number  ( ^ )  of  the  underlying 

3 

lattice.  However  n  could  also  be  larger  than  ^  (See:  Andersen  and  Wheeler  ).  In 
both  cases  the  coordination  number  was  six  (simple  cubic  lattice).  We  will  restrict 
ourselves  also  to  the  simple  cubic  lattice. 

Now  one  can  imagine  two  kinds  of  possible  solutions: 

a)  The  concentrations  of  the  A  molecules  in  the  n  different  ("absolute") 
directions  are  all  equal  and  similarly  for  the  B  molecules  or 

b)  The  concentrations  of  the  A  molecules  in  the  different  absolute 
orientations  are  not  equal  (similarly  for  the  B  molecules),  i.e.  there  is  a  net 
orientation  of  the  molecules  (  liquid  crystals). 

For  the  time  being  we  will  only  consider  the  first  case  (isotropic  solutions)  in 
the  pair  approximation. 

Basic  equations 

The  point  probabilities  (i=l . n)  are  all  equal  since  we  are  dealing 

with  the  isotropic  solution. 


Thus: 


ya-  x» 
*  ”  pT 


where  x^  is  the  concentration  of  the  A  moI-?cules.  A  similar  expression  holds  for  the 


B  molecules,  with  the  normalization: 


(Y  -  Y)= 


Tbe  pair  probabilities  are  given  by  (  in  the  isotropic  solution): 


Variable  (yi)  t 


)  Weight  factor  (g  )  f  Energy  (E. ) 


1  t 


any  pair  of  A  n * 

molecules 

No  arrows  pointing  2(n-l): 

towards  the  other 

molecules 

A  not  pointing  towards  2(n-l) 
the  B  molecule 

B  not  pointing  2(n-l) 

towards  the  A 

molecule 

2 

any  pair  of  B  n2 

molecules 


Normalization: 


Energies  can  be  represented  by  the  matrix  '  *  ;  for  example  L  ~  means  that  U 

3  ^  X  4 

different  from  U2  =  u3  =  U4. 


Relations  between  point  and  pair  probabilities  are: 


*6 +  ^  % 
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It  is  clear  that  here  all  different  are  equal  to  each  other.  However  the 
ratio  of  the  special  contact  point  to  the  sum  of  special  and  normal  contact  points 
for  each  type  of  molecule  (which  should  be  less  or  equal  to  1/*^  depending  on  the 
value  of  n)  might  not  be  conserved.  This  will  be  discussed  later  on. 

(j)  (=^ f  ;  where  f  is  the  free  energy  per  molecule  and  fi  =  1/kT  )  is  then 


given  by 


2. 


§  -  I  51  p  €c  ^  +  xb  ^ ] 


■h 


+ 


is  A0-r^) 

-i-i 

where  \  the  Lagrange  multiplier  for  the  normalization  of  the  probabilities. 
Minimizing  <£  with  respect  to  y^  leads  to: 

y2= 

Now  replacing: 

+  +  y7  =  ^  yi' 

and  putting: 

y2'  =  «vXp(  f  (&A)  (*A  exF  (-  k  ) 

where  K  represents  the  effective  interaction  energy  between  A  and  B  molecules  times 
.  It  follows  immediately  that  K  has  to  satisfy: 

yO’^T*  =■  (n-^ e.  +  (n-i)vT ^UV  ^\  — 7 eT P ^ ^ -t- 

It  can  be  shown  that  this  system  is  now  entirely  analogous  to  a  magnetic  spin-  1/2 
system  with  K/2  replacing  Phase  separation  is  determined  by  the  following 

equation  (for  0  4  0)  : 


k  Si-rvK  (_■$■ 
£  - 


where 


0 


is  given  by: 


The  critical  points  are  determined  by: 
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Kc= 

Phase  separation  occurs  when  K  >  Kc«  Figures  1  to  4  give  the  phase  diagrams  for 
certain  specific  cases.  Since  the  phase  diagrams  are  symmetric  around  x.=.5  only  one 
half  of  the  diagrams  is  plotted.  Also  plotted  in  these  figures  are  the  spinodal 
lines  (denoted  by  s). 

Discussion 

Consider  the  original  Barker  and  Fock  solution.  They  put  constraints  on  the 
different  pair  probabilities  such  that  ratio  of  special  to  all  contact  points  is 
conserved . 


In  our  case  this  ratio  is  not  necessarily  conserved.  To  see  whether  this  ratio 
deviates  much  from  the  exact  ratio,  a  plot  is  made  (fig.  5)  of  this  ratio  (R  divided 
by  the  exact  ratio,  R*=l/6)  vs  temperature  for  two  cases  (corresponding  to  fig.  1 
and'  fig. 2).  The  maximum  deviation  is  seen  to  be  approximately  20%,  which  is, 
considering  that  only  the  pair  approximation  is  used,  not  so  bad. 


case: 


However  our  solution  does  not  give  "unphysical"  phase  diagrams  like  in  their 
for  which  a  phase  diagram  is  given  in  fig.  6. 


i  a 
Z3 


Furthermore,  the  isotropic  solution  also  has  closed  loop  phase  diagrams  for  the  fee 
lattice  i.e.y*12,  whereas  in  Barker  and  F<xk  the  closed  loop  phase  diagrams  are  only 
obtained  for  y  given  by: 


^  -i.vr<^<  4 1-  i\ii 


i 
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in  the  case  of 

In  comparing  the  isotropic  model  with  Andersen  and  Wheeler’s  model  three  main 
differences  can  be  observed. 

I.  In  their  decorated  lattice  model  the  decorated  (or  secondary)  lattice 
points  are  also  occupied  by  molecules. 

II.  In  their  model  only  the  orientation  of  the  molecules  on  the  decorated 
lattice  is  of  importance  (i.e.  no  orientation  of  molecules  on  the 

primary  lattice.) 

III.  Because  of  these  two  features  of  their  model  they  are  able  to  map 
their  model  on  the  spin-1/2  Ising  model,  for  which  the  "exact"  solution  is  known.  It 
can  be  shown  that  in  the  cluster  variation  method  an  increase  in  cluster  of  the  phase 
diagrams.  This  is  illustrated  in  figure  7,  where  the  phase  diagrams  of  Ising  systems 
with  only  repulsive  energies  between  unlike  spins  are  plotted  versus  the  reduced 
temperature  in  the  first  three  approximations. 

It  is  highly  probable  that  in  the  decorated  model,  this  widening  has  also 
occured  because  the  model  is  mapped  on  the  Ising  model  for  which  they  use  the  exact 
solution  (and  not  a  lower  order  approximation). 


Figure  Captions 

l  2 

Fig.  1  Phase  diagram  with  n=6  and  energies  given  by  %  l  »u1=1>  U2=-.25  and  U3=0 

,U=1  and  U2=-.25 


Fig.  3  Phase  diagram  with  n=5001  and  energies  given  by 
Fig.  4  Phase  diagram  with  n=5001  and  energies  given  by 
In  the  above  figures  s  denotes  the  spinodal. 

Fig.  5a  Ratio  of  special  contact  point  to  all  contact  points  of  A  molecules 
normalized  by  the  exact  value  (R*=  1/6)  corresponding  to  fig.  1. 

Fig.  5b  Idem,  but  corresponding  to  fig.  2. 

Fig.  6  Phase  diagram  with  n=6  and  energies  given  by 
(Barker  and  Fock  solution) 

Fig.  7  Phase  diagrams  in  different  approximations 

1.  Mean  field  (Tc=3) 

2.  pair  approximation  (Tc=2.4663) 

3.  4  -  cluster  approximation  (Tc=2.3049) 


U^l, 


u2=.° 


and  U3=-.25 


^11,3^1,  U2=-l  and  U3=-10 
'x\  ,  U1=l  and  U2=-l 


Fig.  2  Phase  diagram  with  n=6  and  energies  given  by  ^ 
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Figure  Captions 

i  2 

Fig.  1  Phase  diagram  with  n=6  and  energies  given  by  2  3  >Uj=l,  U2=-.25  and  U3=0 

,U1=1  and  U2=-.25 


Fig.  3  Phase  diagram  with  n=5001  and  energies  given  by 
Fig.  4  Phase  diagram  with  n=5001  and  energies  given  by 
In  the  above  figures  s  denotes  the  spinodal. 

Fig.  5a  Ratio  of  special  contact  point  to  all  contact  points  of  A  molecules 
normalized  by  the  exact  value  (R*=  1/6)  corresponding  to  fig.  1. 


yil.U^l,  U2=-l  and  U3=-10 
x\  ’  anc* 


Fig.  2  Phase  diagram  with  n=6  and  energies  given  by  ^ 


Fig.  5b  Idem,  but  corresponding  to  fig.  2. 

Fig.  6  Phase  diagram  with  n=6  and  energies  given  by  ^  ,  1^=1,  U?=.0  and  U3=-.25 

(Barker  and  Fock  solution) 


Fig.  7  Phase  diagrams  in  different  approximations 


1.  Mean  field  (T  =3) 

c 

2.  pair  approximation  (T  =2.4663) 

c 

3.  4  -  cluster  approximation  (Tc=2.3049) 
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